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ABSTRACT 

The statistics of gravitationally lensed arcs, which can be used for a variety of cosmological tests, 
are sensitive to the intrinsic shapes of the source galaxies. I present an analytic formalism that makes 
it simple to include elliptical sources in analytic calculations of lens statistics. For cuspy lens models, 
sources with an axis ratio of 2:1 enhance the total number of arcs longer than 10:1 by a factor of order 
two, while modestly decreasing the number ratio of radial arcs to tangential arcs. Source ellipticity 
is therefore an important systematic effect in detailed quantitative studies, but it should not hinder 
cosmological applications such as attempts to constrain cluster dark matter profiles with arc statistics. 

Subject headings: galaxies: clusters: general — dark matter — gravitational lensing 
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1. INTRODUCTION 

A cluster of galaxies can distort background galaxies 
into long, thin arcs by gravitational lensing. The statis- 
tics of lensed arcs can be used to constrain cluster mass 
distributions (e.g., Wu & Hammer 1993; Miralda-Escude 
1993a; Bartelmann, Steinmetz & Weiss 1995; Molikawa 
& Hattori 2000), the population of faint, distant galaxies 
(e.g., Miralda-Escude 1993b; Hamana & Futamase 1997; 
Bczccourt, Pello & Soucail 1998), and cosmological pa- 
rameters (e.g., Wu & Mao 1996; Bartelmann et al. 1998). 
Comparing the number of arcs extended radially relative 
to the cluster center with the number of arcs stretched 
tangentially may be an especially good probe of the cen- 
tral density profile in clusters (Molikawa & Hattori 2000; 
Oguri, Taruya & Suto 2001). 

Most calculations of arc statistics assume that the source 
galaxies are circularly symmetric on the sky. The excep- 
tion is a series of papers by Bartelmann et al. (1995, 1998; 
also Meneghetti et al. 2000), which show that elliptical 
sources increase the likelihood of long, thin arcs but do not 
carefully quantify the effect. As applications of arc statis- 
tics become increasingly detailed, it is important to revisit 
the question of how elliptical sources affect the statistics of 
both tangential and radial arcs. In this paper I develop an 
analytic formalism for including elliptical sources in stud- 
ies of arc statistics, which makes elliptical sources as easy 
to work with as circular sources. Section 2 presents the 
formalism, Section 3 gives examples to quantify the effects 
of elliptical sources, and Section 4 offers conclusions. 

2. ANALYTIC FORMALISM 

The lens equation relates the angular position x of a 
lensed image to the angular position u of the intrinsic 
source, 

u = x V0(x) , (1) 

where </> is the projected gravitational potential of the 
lens, which is related to the projected surface mass density 
£(x) via the two-dimensional Poisson equation V 2 (/>(x) = 



£(x)/£ cr , where S cr is the critical surface density for lens- 
ing (e.g., Schneider, Ehlers & Falco 1992). The distortions 
of an image are determined by the magnification tensor 
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where (f>ij = d 2 4>/dxidxj are second derivatives of the po- 
tential. 

Consider a source with elliptical symmetry, so it is de- 
scribed by a surface brightness distribution with the form 
S(u) = f(p s ), where p s is an ellipse coordinate given by 



where 
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where q s > 1 is the axis ratio of the source ellipse, 9 is 
the orientation angle of the ellipse, R is a rotation ma- 
trix, and t denotes the matrix transpose. Without loss 
of generality, we can choose coordinates aligned with the 
eigenvectors of the lensing magnification tensor, so 9 is the 
angle between the source ellipse and the dominant eigen- 
vector of the magnification tensor (the shear vector). In 
this coordinate system, the magnification tensor can be 
written as 

' qi o 
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where qi > 1 is the ratio of the magnification eigenvalues; 
it gives the axis ratio of an image from a circular source, 
so I refer to it as the lensing axis ratio. Also, a is a scale 
factor that gives the uniform expansion or contraction of 
the image relative to the source. 

Analytic studies of arc statistics often assume a small 
source so the shape of the image is easy to calculate, and 
the assumption appears to be broadly valid (see Hattori, 
Watanabe & Yamashita 1997). In this limit, the image is 
an ellipse described by the surface brightness distribution 



/(x) = S(^i nv -x) = f(pi), where /j,i„ v is the inverse of the 
magnification tensor, and pi is an ellipse coordinate given 
by 



where M 
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(Compare eq. 3.) The axis ratio of the observed image, 
qobs > 1, is given by the eigenvalues of the tensor M to be 



qobs 
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where T = q 2 + q 2 + (qf - l)(q 2 s - 1) COS 2 6 , 
D = q?q 2 s . 



Some simple limits are as follows: 



qobs 








if 6 = tt/2 
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When computing the cross section for arcs, each source 
position must be weighted by the probability that a source 
at that position produces a lensed arc. An arc is defined 
as an image with an axis ratio larger than some thresh- 
old, i.e., q b s > Q for some value of Q; a popular choice is 
Q = 10 (e.g., Wu & Hammer 1993), although Q = A may 
be more relevant for radial arcs (see Oguri et al. 2001). 
Consider a source with axis ratio q s and orientation angle 
9, at a position such that the lensing axis ratio is qi. The 
source produces an arc if 

cos 2 6>Z(Q, qi ,q s ), (13) 

(Q 2 + l)qiq s -Q(qf + q 2 ) 



where £{Q,qi,q. 



Q(q 2 -l)(q 2 -l) 



.(14) 



If the sources have random orientations, the probability of 
obtaining an arc from a source at that position is then 



Parc(Q,qi,Qs 



(2/tt) cos l y/£,(Q,qi,q s 

1 





if < £ < 1 
if e < 
if£>l 

(15) 
Figure 1 shows how to interpret this result. There are 
three simple regions in the (qi,q s ) plane. First, for qi > 
Q q s , the lensing distortion is strong enough that any 
source at that position yields an arc (p arc = 1), regard- 
less of the source orientation. Second, for q s > Qqi, the 
sources are intrinsically long enough that any source at 
that position produces an arc (again p arc = 1), although 
this regime is uninteresting because such elongated sources 
are rare, and because this case is not really lensing. Third, 
for qi q s < Q no source at that position can produce an arc 
(Pare = 0), because neither the intrinsic source shape nor 
the lensing distortion is large enough (for any orientation) . 
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Fig. 1. — The arc probability, p a rc, as a function of qi and 
q s , for an arc threshold Q = 10. Contours arc drawn at p a rc = 
0.95, 0.90, 0.85, . . . starting from the upper right, as indicated. The 
heavy dashed lines indicate p a rc = 1, while the shaded region in the 
lower left indicates p ar c = 0. Source axis ratios larger than q s > 2, 
while implausible, are included to show that p ar c is symmetric in qi 
and q s - 

The remainder of the (qi, q s ) plane has intermediate val- 
ues of Pare, with two interesting qualitative results. First, 
with elliptical sources there are positions where the lensing 
distortion is moderate (qi < Q), but some sources can still 
produce arcs; the arcs are produced when the sources are 
roughly aligned with the lensing distortion. This result 
reproduces the result from Bartelmann et al. (1995) that 
elliptical sources increase the likelihood of long/thin arcs. 
Second, there are positions where the lensing distortion is 
strong (qi > Q), but not all sources yield arcs; sources that 
are approximately orthogonal to the lensing distortion fail 
to produce arcs. 

This analysis makes it simple to include the effects of 
elliptical sources in calculations of arc statistics. Given 
sources with a fixed axis ratio q s , the cross section for arcs 
with threshold Q is 



v{Q,q s ) = dup a 



,(Q,qi(u),q s 



(16) 



where qi (u) is the lensing distortion associated with source 
position u, and the integral extends over the source plane. 
While circular sources yield p arc = or 1 , elliptical sources 
introduce variable weighting into the integral. Because 
eqs. (14) and (15) make the effect easy to compute, ellip- 
tical sources can be added to calculations of lens statistics 
with no more work than circular sources. For source po- 
sitions with multiple images and hence multiple values of 
Pare (one for each image), the largest value of p ar c that 
corresponds to the desired type of arc (radial or tangen- 
tial) should be used. The formalism can be extended to a 
distribution of source shapes. Given a probability distri- 
bution P(q s ) for source shapes, the arc cross section is 

a{Q) = I dn I dq s p ar c(Q,qi(u),q s ) P(q s ) . (17) 



For each source position u, the integral over q s can be 
computed with fast numerical integration. 

3. EXAMPLES 

In proposing that arc statistics can be used to constrain 
the inner profiles of dark matter halos, Molikawa & Hat- 
tori (2000) and Oguri et al. (2001) consider only circular 
sources. I reexamine their proposal to give an example 
of the quantitative effects of source ellipticity. Following 
both studies, consider a dark matter halo with a general- 
ized NFW profile, 
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where ro is a scale radius and po is a characteristic den- 
sity. The density profile has a central cusp p ex r~ 7 ; 
when 7 = 1 it corresponds to the NFW profile (Navarro, 
Frenk & White 1996, 1997), and when 7 = 1.5 it re- 
sembles the steeper cusp advocated by, e.g., Fukushige & 
Makino (1997), Moore et al. (1998, 1999) and Klypin ct al. 
(2000). It is convenient to define a "concentration" param- 
eter C = r2oo/ r *7 where r2oo is the radius within which the 
mean density is 200 times the critical density of the uni- 
verse, often taken to define the edge of halo (e.g., Crone, 
Evrard & Richstone 1994; Cole & Lacey 1996; Navarro et 
al. 1996, 1997); and r* is the radius at which the logarith- 
mic slope of the density is d(\ogp)/d(\ogr) = —2. This 
definition is equivalent to the standard concentration pa- 
rameter for NFW halos, and it is a good generalization 
to 7 7^ 1 (see Keeton & Madau 2001). Like Molikawa & 
Hattori (2000) and Oguri ct al. (2001), I consider spherical 
halos for simplicity, but the formalism presented in Section 
2 is fully general. 

Lcnsing depends on the projected surface mass density 
in units of the critical surface density for lensing. For gen- 
eralized NFW halos, the lensing strength is determined by 
the dimcnsionlcss parameter 

kq = pm = lo4/3 (gW! D^ (1Q) 



(3 



D os 

7)(r2ooAo) 7 ~ 



2-Fi (3-7,3-7,4-7, -r 2 ooAo) 
where M is the halo mass, H is the Hubble parameter 
evaluated at the lens redshift, and 2-Fi(a, b, c, x) is the hy- 
pcrgeometric function (e.g., Press et al. 1992). Also, D i, 
D os , and Di s are angular diameter distances from the ob- 
server to the lens, from the observer to the source, and 
from the lens to the source, respectively. The complete 
equations for generalized NFW lenses are given by, e.g., 
Molikawa & Hattori (2000) and Keeton & Madau (2001). 
As a fiducial model, consider a lens of mass 10 15 h~ l M Q 
at redshift z; = 0.3, with sources at redshift z s = 1, in a 
cosmology with matter density J7m = 0.3 and cosmological 
constant 17a = 0.7. Such halos have a range of concentra- 
tion parameters with median C ~ 4.3 (e.g., Eke, Navarro 
& Steinmetz 2000). The lensing strength parameter is 

then: 

r 0.1776 if 7 = 1.0 

10.0543 if 7 = 1.5 

Increasing or decreasing the strength parameter by 0.2 dex 
corresponds to changes in the halo mass of 0.6 dex (a fac- 
tor of 4), or to changes in the concentration of about 0.18 



Kq = 



(20) 



dex (the la scatter in the concentration; Bullock et al. 
2001). 

Figure 2 shows how arc cross sections depend on the 
source ellipticity e s = 1 — b/a, where a and b are the semi- 
major and semi-minor axes of the source, respectively. In- 
creasing the source ellipticity increases the cross section 
for both tangential and radial arcs, but it affects tangen- 
tial arcs more strongly. From e s = to e s = 0.5 the 
cross section for tangential arcs increases by a factor of 
2.6 for 7 = 1.0 models and 2.3 for 7 = 1.5 models, while 
the cross section for radial arcs increases by 1.7 for both 
models. The enhancement factors are not very sensitive 
to the lensing strength kq. They are somewhat more sen- 
sitive to assumptions about magnification bias (which is 
not included in Figure 2), but factors around 2 are still 
typical. 

The statistic advocated by Molikawa & Hattori (2000) 
is the number ratio of radial and tangential arcs, which is 
shown in Figure 3. Although both types of arcs are en- 
hanced by elliptical sources, radial arcs are affected less 
strongly; hence, the ratio of radial to tangential arcs sys- 
tematically decreases with source ellipticity. From e s = 
to e s = 0.5 the ratio decreases by a factor of about 1.5 
for 7 = 1.0 models and 1.4 for 7 = 1.5 models. The de- 
pendence on the strength parameter kq is weak (which is 
why the number ratio is an attractive statistic in the first 
place). Magnification bias is an important systematic ef- 
fect in the dependence of the arc number ratio on the cusp 
slope 7 (see Oguri et al. 2001), but it does not dramatically 
change the dependence on source ellipticity. 
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Fig. 2. — Arc cross sections versus source ellipticity, computed 
for the fiducial model with a length/width threshold Q = 10. The 
light and heavy curves denote models with 7 = 1.0 and 7 = 1.5, 
respectively. The solid and dashed curves indicate cross sections 
for tangential and radial arcs, respectively. The cross sections are 
computed in angular units and are normalized by the radius r* at 
which the logarithmic slope of the density is —2. 
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Fig. 3. — The number ratio of radial and tangential arcs as a func- 
tion of the source ellipticity, computed with a length/width thresh- 
old Q = 10. The light and heavy curves denote models with 7 = 1.0 
and 7 = 1.5, respectively. The solid curves indicate the fiducial 
model, while the dashed (dotted) curves indicate models where the 
lensing strength is increased (decreased) by 0.2 dcx. 



For modest source ellipticities, e s < 0.5, the change in 
the arc number ratio due to source ellipticity is substan- 
tially smaller than the difference between 7 = 1.0 and 
7 = 1.5 models. Hence, even if the source shape dis- 
tribution is poorly known, it should still be possible for 
arc statistics to distinguish between models with differ- 
ent density cusps. Source ellipticity is thus an important 
systematic effect that must be included to obtain robust 
quantitative results, but it is weaker than effects that one 
may wish to probe with arc statistics. 

These examples are based on axially symmetric lens 
models, as were the studies by Molikawa & Hattori (2000) 



and Oguri et al. (2001). Departures from axial symme- 
try, such as an aspherical cluster potential or subclumps 
in the cluster mass distribution, can affect the number of 
arcs by an order of magnitude or more (Bartelmann et 
al. 1995). While asymmetry in the lens thus appears to 
be more important than non-circular sources, the quanti- 
tative effects of lens asymmetry on the arc number ratio, 
and the combined effects of lens and source asymmetry, are 
under investigation (M. Mencghetti et al., in preparation). 

4. CONCLUSIONS 

I have presented a simple formalism for including the 
effects of elliptical sources in analytic calculations of the 
statistics of lensed arcs. It is straightforward to compute, 
as a function of source position, the probability that a 
source with a given axis ratio and a random orientation 
yields a lensed arc. This probability serves as a weighting 
factor in the integral over source positions to determine the 
lensing cross section. The formalism is fully analytic, so el- 
liptical sources require no more work than circular sources; 
and it is easily extended to a distribution of source shapes. 

In generalized NFW lenses, elliptical sources enhance 
the likelihood of producing both tangential and radial arcs; 
sources with ellipticity e s = 0.5 increase the expected num- 
ber of arcs by about a factor of two (relative to circular 
sources). Radial and tangential arcs are affected differ- 
ently, so the number ratio of radial arcs to tangential arcs 
decreases as e s increases; from e s = to e s = 0.5 the ratio 
drops by a factor of about 1.5. This change is an important 
systematic effect in quantitative studies of arc statistics, 
but by itself it should not significantly limit attempts to 
constrain the dark matter profile of clusters. These con- 
clusions are drawn from examples with spherical lenses, 
and it will be interesting to see whether elliptical sources 
have any different effects when the lenses are non-spherical 
(M. Meneghetti et al., in preparation). 
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